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ABSTRACT 

The width (Ai^), root mean squared amplitude (rg) of lower and upper kHz quasi-periodic 
oscillations (QPOs) from accreting neutron stars vary with frequency. Similarly, the QPO fre- 
quency varies with the source count rate (S). Hence, the significance of a QPO, scaling as 
S X r^/ V Aiy will also depend on frequency. In addition, the significance also scales up with 
the square root of the integration time of the Fourier power density spectrum (T). In prac- 
tice and in most analysis, T is constrained (e.g. limited by the RXTE orbital period if the 
source is occulted by the earth) or minimized to increase the number of detections. Conse- 
quently, depending on the way data are considered, kHz QPOs (the lower, the upper or both 
kHz QPOs) may be detected only over a limited range of their frequency spans or detected 
predominantly at some frequencies, leading potentially to biases in the observed distributions 
of frequencies or frequency ratios. Although subject of much controversy, an observed clus- 
tering of QPO frequency ratios around 3/2 in Sco X-1, also seen in other sources, has been 
previously used as an argument supporting resonance based models of neutron star QPOs. 
In this paper, we measure how the statistical significance of both kHz QPOs vary with fre- 
quency. For this purpose, we consider three prototype neutron star kHz QPO sources, namely 
4U1636-536, 4U0614H-091 and Sco X-1, whose QPO parameters have different, though rep- 
resentative, dependence with frequency. As the significance of QPO detection depends on 
frequency, we show that in sensitivity-limited observations (as in the case of the RXTE/PCA), 
a simultaneous detection of both the lower and upper kHz QPOs can only be achieved over 
limited frequency ranges. As a result, even a uniform distribution of QPO frequencies will 
lead to peaks (in particular around 3/2) in the histogram of ratios of simultaneous kHz QPO 
frequencies. This implies that the observed clustering of ratios of twin QPO frequencies does 
not provide any evidence for intrinsically preferred frequency ratios in those systems, thus 
weakening the case for a resonance mechanism at the origin of neutron star kHz QPOs. 
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1 INTRODUCTION 

Despite a wealth of observational data, accumulated since their dis- 
covery in 1995, there is not yet a commonly accepted model of 
neutron star kHz QPOs. In recent years, a lot of efforts has been 
devoted to studying the distribution of the observed frequencies (or 
frequency ratios), in an attempt to provide insights to theoretical 
models. This was particularly relevant for resonance based mod- 
els which predicted that both the ratio of frequencies as well as 
the frequencies themselves should cluster around specific values. 
Apparent support for these models came from the observation that 
the ratios of simultaneous lower and upper QPO frequencies from 
Sco X-1 clustered around 1.5, su ggestive of a 3:2 resonance mecha- 
nism jAbramowicz et al .l2003al) . Additional support came simulta- 
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neously from the observations of bl ack hole high fre quency QPOs 
associated with fixed 3/2 ratios (e.g. lRemillardll2005l . for a review). 
The clustering of frequency ratios, claimed in other neutron star 
systems as well (Toroketal. 2 008al . for 4U1636-536), remains a 
controversial issue sinc e then. In particula r , the si gnificance of the 
cl ustering presented by Abramowicz et al.l ( l2003al) was challenged 
bv lBelloni et"ZI ilOOsS . 

The signal to noise ratio rio- at which aperiodic variabil- 
ity such as a kHz QPO is detected in a photon counting experi- 
ment is approxim ately ( iLewin. van Paradiis. & van der Klis|[T983 : 
Ivan derKlisI 19891) : 
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where S and B are source and background count rate, respectively, 
rs is the root mean squared amplitude of the variability expressed 
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as a fraction of S, T the integration time and Av the. bandwidtfi 
of tiie variability. Hereafter, we will neglect B over S, because 
for most QPO sources, B (a few te ns of co unts/s/PC U, where 
PCU s tands for PCA Unit, 5 at most, iBradt. Rot hschil d & SwanlJ 
( Il993h ) is at least a factor of 10 lower than 5", which is typically 
several hundreds of counts/s/PCU. In equationfl] rg and Av vary 
with frequency (ly) tearret. Olive. & Milled [20061) . hence ricr will 
depend also on v. The frequency dependence of rg and Ai/ is a 
non monotonic function of frequency and has a different shape for 
the lower and upper kHz QPOs. In addition, in a given system, on 
average larger luminosities (hence count rates, S) correspond to 
larger v jpord e t al. 2000). Actually, the situation is more complex, 
because of the parallel track phenomenon, which shows that i/ is 
not uniquely determined by S, but that the spr ead of S over whic h 
QPOs are detected increases wit h frequencv |van der Klisll200ll) . 
This is illustrated in Figure 2 in iBarret. Olive. & Millej hoOSst) 
where, in a frequency versus count rate diagram, it is mostly the 
left upper part of the diagram that is populated. The implications 
of the above are twofold. First, the integration time at which a kHz 
QPO is detected above a certain significance threshold will depend 
on frequency. Second, if T is constrained, kHz QPOs, either the 
lower, the upper or both may reach the significance threshold only 
over limited frequency ranges. This directly implies that the ob- 
served distribution of frequencies or frequency ratios will strongly 
depend on the way the data are considered. 

When dealing with frequency distributions, two main types 
of analysis have been applied so far, both suffering from intrin- 
sic biases. Methods consisting of averaging large amount of data, 
based on some selection crit eria, e.g. spectral colors, position along 
a color-color diagram (e.g. Ijonker et al. 2002, for GX5-1), same 
type of QPO falling in a given frequency interval, as for the shift 
and add procedure (e.g. iMendez & van der Klisll2OO0l . for Sco X- 
1), erase the underlying distribution of frequencies and cannot be 
considered for this purpose. The first applicable method consid- 
ers observations of similar durations. This is the case in particular 
when data are analyzed ObsID per ObsID, whose durations are lim- 
ited by the RXTE orbital period to typically ~ 2 — 3 kseconds (by 
definition, an ObsID corresponds to a single observation, where ob- 
servation refers to a temporally contiguous collection of data from a 
single pointing). This method will obviously bias the distribution of 
frequencies or ratios at frequencies where the two QPOs can be de- 
tected simultaneously on a typical ObsID timescale. Such a m ethod 
was recently applied to 4U1636-536 by iTorok et alj l l2008al) who 
reporte d a clustering of ratios at 3/2 and possibly at 5/4, to 4U1820- 
303 by IB arret & Boutelieil ( 1200 8l) who found a clustering of ra- 
tio s of simultaneous twin QPOs aro und 4/3 and to 4U06 14-1-091 
bv lBoutelier. Barret. & Miller! ( |2009|) . In both studies, it was found 
that the distribution of ratios computed from the single lower QPO 
frequencies, estimating the upper QPO frequencies from their lin- 
ear relationship differed from the distribution of ratios of simul- 
taneous twin QPOs, as expected if the condition for detecting the 
lower QPO differed from the one for detecting the two QPOs si- 
multaneously. 

As to increase the number of detected frequencies, the sec- 
ond method minimizes the PDS integration time to detect one of 
the two QPOs (generally the lower QPO) and computes the miss- 
ing frequency through the linear relationship linking the lower and 
upper QPO frequencies. This method will bias the distributions at 
frequencies where the lower QPO is easy to detect (e.g. when its 
R MS amplitude / quality factor reaches a maximum). It was used 
bv lBelloniet"ai] ( l2005bl) who averaged a variable number of PDS 
(with integration time ranging from 64 seconds to 1280 seconds) 



to obtain a significant detection of the lower QPO. They found that 
the distribution of lower QPO frequencies showed a pronounced 
peak around 850 Hz in 4U1636-536. Other less prominent peaks 
were present in the di stribution of f r equen cies, all centered on dif- 
ferent ratios. Later on. lBelloni et aLlfcoOVl) monitored 4U1636-536 
for 1.5 years, during regular 2 kilosecond pointings separated by 
2 days, and found that the distribution of lower QPO frequencies 
differed from the previous one. In particular, the frequency distri- 
bution of the lower QPO extended down to 600 Hz and the overall 
distribution was flatter than previously measured. They also noticed 
that when considering one frequency per segment of 2 kseconds, 
the upper kHz QPO was detected only between 520 Hz and 950 Hz 
with an even distribution. All this led lBelloni et alj | |2007|) to con- 
clude that the kHz QPOs in 4U1636-536 did not show intrinsically 
preferred frequencies. 

In this paper, we measure the frequency dependence of and 
show how it affects the distribution of frequencies measured with 
the two methods discussed a bove. This paper extends on previous 
work bv lT6r6keta"i] ( l2008d) who did not directly compute Ua- but 
rather scaled its value as the product of rg/V Au: the frequency 
dependence of rs and Au was recovered for 4U1636-536 by inter- 
polating the values reported by[Barret, Olive, & Miller (2005i), i.e. 
derived from the analysis of a subset of the data considered here. 
Here, we have chosen three systems 4U1636-536, 4U0614-I-091 
and Sco X-1, as being representative of low and high luminos- 
ity neutron star QPO sources, both in terms of RMS amplitudes 
and quality factors. Both sources have been extensively observed 
by RXTE; several hundreds of ObsIDs spanning over more than 
10 years, thus minimizing any biases in the frequency distribu- 
tion of their QPOs that could arise as a result of a non uniform 
sampling of the source states. This matters because iBelloni et al.l 
12005b.) have shown that a time limited and sparse sampling of 
a source, whose QPO frequency was modeled by a random walk 
could pro duce peaked frequen cy distributions. In addition, as dis- 
cussed bv lBelloni eTaklfeOOTh . in the case of 4U1636-536, the long 
term evolution of the source intensity may also affect the observed 
frequency distributions. 

In the following section, we present the histogram of fre- 
quency ratios, as derived from the homogenous processing of 
all archival data available to date for both 4U1636-536 and 
4U06 1 4-1-09 1 . We then describe the analysis scheme used to recover 
over the full frequency span of both the lower and upper QPOs. 
Then, knowing ria as a function of frequency, assuming an uniform 
distribution of frequencies for the upper QPO, we show that we can 
reproduce, without any further assumptions, the clustering of fre- 
quency ratios seen in the data. We also illustrate how the method 
of optimizing the PDS integration time introduce strong biases in 
the histogram of ratios. Finally, we then extend our analysis to Sco 
X- 1 before concluding. 



2 DATA ANALYSIS 

We have first considered all archival data from the Rossi X-ray 
Timing Explorer for 4U1636-5 36 and 4U0614-I-091. Our anal y 



U ur anal y- 
sis bu i lds up on previous work bv lBarret. Olive. & Mille3 (l2005j|bl 
l2006l |2007|) for 4U1636-536 and bv iBoutelier Barret. & MilleJ 



20091) for 4U0614-I-091. Data are considered per ObsID. For each 



ObsID, we have computed 16 second PDS with a 1 Hz resolu- 
tion, using events record ed between 2 and 4 keV. The PDS are 
normalized according to iLeahv et al.l l ll983l) , so that the Poisson 
noise level is constant around 2. As previously shown for 4U1636- 
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Figure 1. Histograms of ratios estimated from tlie frequencies of all single lower QPOs, all lower in twin QPOs (area in grey), all single upper QPOs, all upper 
in twin QPOs (ai'ea in grey) and all twin QPOs. In the case of a single QPO, the frequency of the other was computed from the linear function linking the two 
frequencies. 4U1636-536 is shown on the left and 4U0614+091 on the right hand side. Each frequency in that plot corresponds roughly to a similar integration 
time of the PDS (a typical ObsID duration). 



536, the lower QPO can be so strong and so narrow that in some 
ObsIDs, the frequency drift within the ObsID can be corrected 
for using a tracking procedur e similar to the one described in 
iBarret. Olive. & Millej f2005ij) . Whenever possible, all the 16 sec- 
ond PDS are then aligned to the mean frequency of the lower QPO 
within the ObsID, to produce one single ObsID averaged PDS (oth- 
erwise they are averaged directly). This enables us to obtain mean- 
ingful parameters (width and amplitude) averaged over the ObsID, 
helping us to detect the upper QPO in some case, and in all cases, 
to identify wether it is a lower or an upper kHz QPO (see below). 
In 4U0614-I-091, the QPO is not strong enough to enable such a 
analysis. The ObsID averaged PDS is then blindly searched for ex- 
cess power between 3 00 Hz and 1400 Hz using a scanning tech- 
nique, as presented in iBoirin et al. I 1 I2OOOI) . Each excess (at most 
the 2 strongest) is then fitted with a Lorentzian with three free pa- 
rameters; frequency, full width at half maximum (constrained to 
range from 2 to 1000 Hz), and amplitude (equal to the integrated 
power of the Lorentzian). The Poisson noise level is fitted sepa- 
rately above 1400 Hz and then frozen when fitting the QPOs (the 
fitted noise level is generally close to 2, indicating no significant 
deadtime). Errors on each parameter are computed with Ax^ ~ 1- 
As in previous papers in this field, our threshold for QPOs is re- 
lated to the ratio (hereafter R) of the Lorentzian amplitude to its la 
enw (R was often quoted and used as a significance). In this paper, 
our threshold for i? is 3, meaning that we consider only QPOs for 
which we can measure the power of the Lorentzian with an accu- 
racy of 3(T or more. Such a threshold corresponds to a ~ 5.5 — Go- 
excess power in the PDS for a single trial, equivalent to ~ 4(7 sig- 
nificance i f we account for t he number of trials of the scanning 
procedure l IvanderKlisI 19891) . 

2.1 Histograms of ratios 

At the end of this first stage of our analysis, we have a set of ObsID 
averaged PDS (i.e. of comparable durations) with one or two sig- 



nificant QPOs in. For those averaged PDS in which only one signif- 
icant QPO is detected, we identify whether it is a lower or an upper 
by placing the QPO parameters in a RMS-:^ and quality factor- 
V diagram following ISarret. Olive. & Millet! 1 I2OO6I) . In addition to 
the histogram of ratios derived from simultaneous twin QPO fre- 
quencies, using the nearly linear relationship between the lower and 
upper QPO frequencies, it is possible to compute the histograms 
of ratios estimated from the frequencies of all single lower QPOs 
and all single upper QPOs. The histograms are presented in Figure 
[TJ the frequencies entering the histograms are the m ean frequen- 
cies ove r the ObsID. The clusters previously found bv lTorok et al.l 
l l2008j) around 1.5 and 1.25 for 4U 1636-536 are reproduced and 
clearly the histograms computed from either the single lower or 
single upper QPOs is different from the one computed from simul- 
taneous twin QPOs. A cluster of simultaneous twin QPOs centered 
around 1.5 is also seen from 4U0614-I-091. 

2.2 Frequency dependence of QPO significance 

We now wish to estimate how the significance of kHz QPOs varies 
with f requency. For this pu rpose, we apply the shift-and-add tech- 
nique ( lMendezetalJl998h . As originally discussed, this technique 
enables us to combine large amount of data, so that both the lower 
and upper QPOs can be detected simultaneously (this is possible 
for some ObsID, but generally requires longer integration times of 
the PDS). We shift-and-add all ObsID averaged PDS, considering 
all those with a QPO of the same type falling in the same frequency 
interval. The width of the frequency interval has been chosen to 
group enough PDS as to detect the two QPOs simultaneously, with 
a high significance. The large amount of data in the RXTE archive 
makes possible to consider interval of 30 Hz for 4U1636-536 and 
50 Hz for 4U06 14-1-091. Obviously, not all frequency intervals cor- 
respond to the same total PDS integration time, so we normalized 
the significance using equation[T]to the square root of the observing 
time of a typical ObsID, which we assumed to be 3 kseconds (this 
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Figure 2. The significance of the PDS excess power in 4U1636-536 (left) and 4U0614+091 (right) against frequency, as recovered after applying the shift- 
and-add technique and norm alizing the significa nce to the same observing time, assumed to be 3 kseconds, i.e. comparable to a typical ObsID exposure. The 
significance is compute from lBoirin et al.1 i200d) . It is equivalent to the significance given in equation 1 with a 1/4 scaling coefficient. The horizontal dashed 
line indicates the 5.5 sigma (single trial) statistical significance threshold used in our analysis. The same way that not all frequency bins correspond to the same 
PDS integration time, not all frequency bins correspond to the same number of active PCUs. To account for this, the significance has also been normalized to 
the mean number of active PCUs over the whole frequency span (typically around 4). Neither the lower or upper QPO can be detected over their full frequency 
span on such a timescale, con sidering a threshold of S.Scr. In particular, the gap in t he distribution of frequencies of the upper QPO visible in 4U1636-536 
iBarret. Olive. & MilleJIlOOd) and in 4U06I4+091 iBouteher, Barret, & Milletll2009h is now explained, as being caused by the significance being lower than 
the threshold generally used in previous analysis. We use the same scale of the y axis to show that the lower QPO is much easier to detect in 4U1636-536 than 
in 4U06 14+091. 



is about the mean duration of all ObsIDs from 4U0614+091: it is 
closer to 2.5 kseconds for 4U1636-536). 

The result of such analysis is shown on Figure |2] for both 
4U1636-536 and 4U0614+091, which shows the frequency depen- 
dence of the QPO significance (excess power in the PDS, single 
trial) normalized to the same integration time. Figure [2] demon- 
strates that neither the lower or upper QPO are detected over their 
full frequency span, for a commonly used single trial significance 
threshold of 5.5a. As discussed earlier, na depends on S, rs and 
Av, how those quantity vary with frequency for the lower kHz 
QPO in 4U1636-536 is shown in Fig.[3] The maximum of corre- 
sponds to the maximum of rs and the minimum of Az^ (or the max- 
imum of the QPO coherence). From the lowest QPO frequency to 
the highest one, the mean count rate increases by less than a factor 
of 2. 

For the upper QPO, as shown in Figure [3] is rather flat: 
the decrease of rs is compensated by a decrease of Ai/ and an 
increase of S. There is a l ocal minimum around 1100 Hz in both 
sources. For 4U0614-F091, lBoutelier. Barret, & Millerl ( l2009h have 
shown that this was related to a local decrease of rs. For 4U 1636- 
536 (and also for 4U1820-303, Ustrret & Boutelie:| ( |2008|) ), it is 
associated with a local increase of the width of the upper QPO 
(see Fig.[3]l. This local minimum provides an explanation for the 
gap in the distr ibution of frequencies of the upper QPO visible 
in 4U1636-536 fearret, Olive, & Mi llerll2006l) and in 4U0614-^091 
teoutelier. Barr et, & Miller 2009) (when data are analyzed per Ob- 
sID). After the local minimum, reaches a local maximum, en- 
abling us to detect twin QPOs over a narrow range of frequencies, 
which in turn produces a second cluster of ratios, e.g. below 1.3 in 
4U1636-536. 



3 SIMULATIONS 

Now we wish to visualize the influence of the frequency depen- 
dence of the significance of the QPO on the observed distribution 
of ratios. For this, we assume an uniform distribution of frequen- 
cies for the upper QPO which spans a wider range of frequencies 
than the lower QPO. For example, for 4U1636-536, the upper QPO 
frequency spans from 550 Hz to 1250 Hz, while the lower QPO 
frequency varies only from 550 Hz to 950 Hz. Obviously, by draw- 
ing a uniform distribution of frequencies within the range of the 
lower QPO frequency, one would miss all the single upper QPO 
frequencies below ~ 800 Hz (those QPOs populate the right part 
of the histograms of ratios predicted from the upper QPO also, mid- 
dle panel in Figure[TJ. For each upper QPO frequency so generated, 
we compute the frequency of the lower QPO, from a linear function 
z^uppcr = a X t'lowcr + b with a = 0.70 and b — 520 for 4U1636- 
536 and a = 1.0, and b = 320 for 4U0614-I-091 respectively. The 
a and b parameters were derived for each source by fitting twin 
QPO frequencies detected simultaneously over the same ObsID. 
We have checked that the parameters are consistent with each other 
when fitting the twin QPO frequencies recovered from the shift- 
and-add procedure described above. Instead of fitting with a linear 
function, we have also considered a power law function, leading 
to very similar results, as expected since the two functions are es- 
sentially identical in the range where the two QPOs can be detected 
simultaneously. We discard those predicted lower QPO frequencies 
falling outside the range of detected frequencies, i.e. below 500 Hz 
and above 950 Hz. We interpolate between the values of Uo- shown 
in Figure [2] to estimate the corresponding significance for all gen- 
erated lower and upper QPO frequencies. 
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Figure 3. The RMS amplitude (rg), the width (Ai/), and the mean source count rate (normalized by the mean number of operating PCUs) for the lower (left) 
and upper (right) QPO in 4U1636-536, as averaged over frequency bins of 50 Hz and 75 Hz respectively. This is the combination of these three parameters as 
described in equation 1 which explains the dependence of with frequency shown in Figure|2] The local minimum of around 1 100 Hz visible in Figure 
|2]for the upper QPO is related to a local increase of the QPO width. 



3.1 Histograms with a constant integration time 

We first consider the case in which the significance of each QPO is 
computed for a similar integration time (3 kseconds), i.e. reproduc- 
ing the case where the real data are analyzed ObsID per ObsID (the 
frequency generated represents the mean QPO frequency within an 
ObsID). The results of this simulation, considering a significance 
threshold of 5.5ct, are shown in Figure|4]in a form similar to Figure 
12] What is striking from this figure is that, even this oversimplified 
simulation, reproduces the shape of the histograms and the cluster- 
ing of ratios seen in the data. For example, the second cluster of 
ratios of twin QPOs in 4U1636-536 below 1.3 corresponds to the 
local maximum of the curve for the upper QPO around 1150 
Hz. Finally, it is obvious that the distributions of ratios computed 
from one of the two QPO frequencies are also not complete and 
not representative of the underlying distribution. The good agree- 
ment between the simulations and the data clearly indicates that the 
observations performed by RXTE over the ~ 10 year period have 
provided a close to uniform sampling of the kHz QPOs over their 
entire frequency span. 



3.2 Histograms with a variable integration time 

We have simulated again a uniform distribution of lower QPO fre- 
quencies, generating 590 random frequencies corresponding to 590 
ObsIDs of 3 kseconds (as is the case for 4U1636-536). Using Fig- 
ure [2] and equation 1, it is possible to compute the minimum PDS 
integration time, required for the QPO to exceed a significance 
threshold of 5.5(7, hence the number of detections associated with 
each segment of 3 kseconds (e.g. for 4U 1636-536, at 850 Hz, the 
significance reaches ~ 40a for 3 kseconds: this corresponds to 
5.5(T in about 60 seconds, hence a segment of 3 kseconds will con- 
tribute 50 frequencies around 850 Hz in the histogram). Having 
the lower QPO frequencies, the upper QPO frequencies are esti- 
mated with the linear function described above. From this, one can 



produce an histogram with all the detections, keeping in mind that 
not all detections correspond to the same PDS integration time. An 
example of such simulation is shown in Fig[5] The histogram of ra- 
tios is strongly peaked around 1 .3. as found in data analyzed by the 
same method by Belloni et al. (2005b): the main peak was found 
at 1.30±0.02 for — 856 ± 2. Such histogram (which just re- 
flects the variation of ric, with frequency) is much more biased than 
the histogram of single lower QPO computed from detections over 
similar integration times (see Figure|4] left hand side, second panel 
from the top). 

3.3 More realistic simulations 

We have also added complexity to our simulations by modeling 
the parallel-track phenomenon (or the frequency dependence with 
count rates), the spread of the ObsID exposure times around their 
mean value and a varying number of active PCA units (the latter 
two quantities being obviously not frequency dependent). We have 
also included a jitter in frequency along the linear correlation be- 
tween the two frequencies. In all cases, this does not change signifi- 
cantly the histograms simulated, it just increases slightly the spread 
of the clusters, making them more similar with the data. This is il- 
lustrated for the case of 4U0614-I-091 in Figure[6] The gap between 

1.4 and 1.5 in the ratio distribution is now filled in. 



4 SCOX-1 

Building on the above findings, it is worth looking back at the 
historical case of Sco X-1, and see whether the same effect can 
provide an explanation fo r the clustering of ratio claimed at 1.5 
dAbramowicz et al ] |2003al) . For this purpose, we have reprocessed 
all the a vailable data from the RXTE archive, extending on t he dat a 
used bv lvan der Klis et all 119971) : iMendez & van der Klii ilOOd) . 
We have used all science event and binned mode data with a time 
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Figure 4. Simulated histograms of ratios expected from a uniform distribution (top panel, assuming that the upper QPO frequency distribution is uniform 
between 550 Hz and 1250 Hz), and the one derived from the frequencies of all single lower QPOs (second panel from top, the shaded area recalls the histogram 
of ratio from simultaneous twin QPO frequencies shown at the bottom), all single upper QPOs (third panel from top, the s haded a rea recalls the histogram 
of ratio from simultaneous twin QPO frequencies), and all twin QPOs (bottom panel). As first pointed out bv lBelloni et"^ i2005bl) , as a consequence of the 
nearly linear relationship between the two QPO frequencies, note that the histogram of ratios from a uniform distribution of frequencies is peaked towards 
lower ratios. In the case of a single QPO, the frequency of the other was computed from the linear function linking the two frequencies. 4U1636-536 and 
4U06 14+091 ai'e shown on the left and right hand side respectively. The main peaks as well as the secondary peaks of the observed distributions are nicely 
reproduced with the simulations, even with our simplified assumptions. 



izr ' ■ 1 p ' ■ ' n F — ' — ' — I — ' — ' — ' — I — ' — ' — ' — r 




QPO Frequency (Hz) 



Figure 7. Left) The histogram of ratios from all lower, all upper and all twin kHz QPO frequencies. In the top two panels, the histogram in grey is the one 
derived from simultaneous twins. As can be seen, twin QPOs are generally detected simultaneously for Sco X-1, and there are no occurrence where the lower 
kHz QPO is detected alone (unlike in the two other sources). In the middle panel, note that there is one single detection of a single upper QPO, at a frequency 
just above ~ 750 Hz, con'esponding to a ratio ~ 1.8 (this is the lowest upper QPO frequency reported for the source so far). The gap in the ratio distribution 
corresponds to a frequency gap about 100 Hz. Right) The significance of the PDS excess power in Sco X-1 against frequency, as recovered after applying 
the shift-and-add technique over frequency intervals of 50 Hz and normalizing the significance to the same observing time, assumed to be 3 kseconds, i.e. 
comparable to a typical ObsID exposure. The lower QPO significance is always lower than the upper QPO significance. The lower QPO is detected over a 
relatively nan'ow range of frequencies about ~ 300 Hz wide. The first bin of the upper QPO histogram, where the significance reaches a maximum, contains 
only one ObsID (the upper QPO frequency is at about 750 Hz) 
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Figure 5. Simulated histogram of ratios estimated from tlie frequencies of 
all single lower QPOs, assuming a uniform distribution of frequencies, and a 
variable integration time of the PDS for detecting the QPO. The histograms 
peaks at the ratio which con'esponds to a frequency of ~ 850 Hz for the 
lower QPO. i.e. where it can be followed on the shorter timescales. The 
shaded region delimited by a dashed line corresponds to the ratio histogram, 
expected from the uniform distribution, and considering one frequency per 
observation. 
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Figure 6. Simulated histogram of ratios for 4U06 14+091 for a more sophis- 
ticated simulation, which accounts for the fact that the spread of count rates 
at which QPOs are detected increases with frequencies (the parallel tracks), 
for a spread in the exposure time between 2 and 4 kseconds, for a spread in 
the number of active PCUs (between 3 and 5) and finally for a 10 Hz jitter 
in the relation linking the two QPO frequencies. As a result, the simulated 
clusters show a wider spread more in line with the real data (see Figure[T]) 



resolution better than 250 /is. We have again computed 16 sec- 
ond PDS with a 1 Hz resolution, using events recorded between 
2 and 40 keV. The analysis here is complicated by deadtime effect, 
which we have to account for in the PDS analysis. We have mod- 
eled the modification of the PDS with a power law of index close 
to zero at high frequencies, and fitted the QPOs with the method 
described above. Over our data set, we have detected about 100 
QPOs from the source on a typical ObsID timescales. Likewise for 
the two other sources, we identify single QPOs by placing them 
on a RMS/quality factor versus frequency plot. The lower QPO 
frequency varies from 550 Hz to 800 Hz, while the upper QPO 
frequency varies between 750 Hz and 1100 Hz. Previo us analysis 
jvan der Klis et alJI 19971 : iMendez & van der Klisll2000l) . based on 
a smaller data set, found the upper QPO frequency to range from 
about 850 Hz to 1150 Hz. We report for the first time, an upper 
QPO at a frequency of 750 Hz, detected in the ObsID 93067-01- 
01-02 recorded on July 3rd, 2007 at 10:33 am. Next to the 750 Hz 
upper QPO, the second lowest upper QPO frequency is at ~ 850 
Hz, implying a gap of 100 Hz in its frequency distribution. The his- 
tograms of ratios are shown in Figure [7] As previously found the 
histograms cluster around 1.5. We have then applied the shift-and- 
add technique to the ObsID averaged PDS as described above, and 
computed the average statistical significance of both QPOs, which 
we then normalized to a PDS integration time of 3 kseconds. The 
results are also shown on Figure |7] In Sco X-1, the lower QPO 
is detected on such a timescale only in a narrow 300 Hz interval, 
whereas the upper QPO is always more significant and detected 
over a slightly wider frequency span ~ 350 Hz wide (note however 
that there is only one ObsID in the first histogram bin where the 
significance reaches a maximum). Note that the narrowness of the 
frequency range for the both QPOs limits by itself the range of pos- 
sible ratios between 1.3 and 1.6. Figure[8]shows how S, rs and 
depend on frequency for both the lower and upper kHz QPOs. As 
can be seen, the 750 Hz QPO follows the same trends rs and Av 
as upper QPOs detected above 800 Hz, and is therefore very likely 
a single upper QPO. Note also that the mean count rate varies by 
less than 20% of the mean, across the frequency range over which 
both QPOs are detected. 

As for the two other sources, in an attempt to compare with 
the observations, we have simulated a uniform distribution of up- 
per QPO frequencies between 850 Hz and 1100 Hz wher e the 
bulk of the upper QPOs are detected jvan der Klis et al ]| 19971) . We 
have computed the histogram of ratios using a linear function with 
a = 0.81 and b = 416, which is equivalent to a power law relation 
i/iower = A(i^uppcr/1000 Hz)^, with A= 717 and B= 1.78 in the 
range of frequency considered. We have also included a spread in 
the count rate, exposure time, number of PCUs and a jitter in fre- 
quency along the correlation line. The results are shown in Figure 
|9] Clearly, the simulations reproduces nicely the observed distribu- 
tions, in particular the strong clustering of ratio around 1.5. We thus 
conclude that as in the two other sources, the apparent clustering of 
frequency ratios in Sco X-1 is also consistent with a uniform dis- 
tribution of frequencies between 850 Hz and 1 100 Hz. However, in 
Sco X-1 the clustering seems to be amplified compared to the other 
sources, by the fact that both QPOs can be detected simultaneously 
over an narrower frequency range. 



5 CONCLUSIONS 

We have shown that, because of the way the significance of kHz 
QPOs evolves with frequency, as a result of the QPO RMS ampli- 
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Figure 8. The RMS amplitude (rg), the width (Ai/), and the mean source count rate (normalized by the mean number of operating PCUs) for the lower (left) 
and upper (right) QPO in Sco X- 1 , as averaged over frequency bins of 25 Hz. This is the combination of these three parameters as described in equation 1 
which explains the dependence of rio- with frequency shown in Figure|2] 



tude and width being frequency dependent quantities, in sensitivity 
limited observations the observed clusters in frequency ratios of 
simultaneous twin QPO frequencies, in particular around 1.5, are 
consistent with an underlying uniform frequency distribution. In 
other words, this implies that the clustering of ratios does not pro- 
vide any evidence for preferred fre quency ratios i i i those systems. A 
similar conclusion was reached bv lBelloni et al.l ( |2007|) who found 
that the distribution of frequencies could be significantly different 
from one data set to the other. The clustering is a consequence of the 
fact, that albeit always present, the lower and upper QPOs can only 
be detected together on the same integration timescales over a lim- 
ited frequency range. The clustering of ratios is more pronounced 
in Sco X-1 because the range of simultaneous detection of the two 
QPOs is narrower (~ 250 - 300 Hz instead of 400 Hz). This con- 
clusion applies to a sample of systems representative of all types 
of QPO sources. As said above, clustering of ratios caused by an 
incomplete sampling of the sourc e states, whether or not its QPO 
frequency follows a random walk jSelloni et alj2005bl) is more and 
more unlikely, as more and more data become available. 

Any models, attempting to reproduce the clusters of ratios, 
must therefore predict first a nearly linear relationship between 
the two QPO frequencies in the range where the bulk of QPOs 
are detected, a s most models do, includin g non-linear resonance 
based models jAbramowicz et all l2003bl) . and second, a set of 
well matched QPO parameters (RMS amplitude and width) which 
allows the two QPOs to be preferentially detected on similar 
timescales over limited frequency ranges. IJorok (2009) already 
pointed out that at frequencies corresponding to the 3/2 ratio, the 
RMS amplitudes of the lower and upper QPOs as measured in the 
2-40 keV band are equal (this is where the significance of the two 
QPOs are comparable, around 600 Hz, 900 Hz for the lower and 
upper respectively, see Figure |2j. In the framework of the non - 
linear resonance model proposed by lAbramowicz et aL I j2003bl) . 
this has recently been discussed as a possible signat ure of a res- 
onant energy exchange between the two QPO modes dHorak et al.l 
I2OO9I) . Although such a possibility should be investigated further, it 



is worth stressing that non linear resonance was invoked for neutron 
star QPOs by interpreting the 3/2 clustering of frequency ratios as 
reflecting a non uniform distribution o f frequencies along the lin e 
of correlation predicted by the model jAbramowicz et al]|2003bl) . 
This is the reason why the results presented in this paper weaken 
the case for a resonance mechanism at the origin of neutron star 
kHz QPOs, as they indicate that even a uniform distribution of fre- 
quencies will produce peaks in the ratio distribution around 3/2. 
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